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We construct models of universe with a generalized equation of state p = (ap + kp 



l+l/n\ 2 



having 



a linear component and a polytropic component. The linear equation of state p = ape 2 with 
— 1 < a < 1 describes radiation (a = 1/3), pressureless matter (a = 0), stiff matter (a — 1), and 
vacuum energy (a = —1). The polytropic equation of state p = kp 1+1 ' n c 2 may be due to Bose- 
Einstein condensates with repulsive (k > 0) or attractive (k < 0) self-interaction, or have another 
origin. In this paper, we consider the case where the density increases as the universe expands. This 
corresponds to a "phantom universe" for which w = p/pc 2 < — 1 (this requires k < 0). We complete 
previous investigations on this problem and analyze in detail the different possibilities. We describe 
the singularities using the classification of [S. Nojiri, S.D. Odintsov, S. Tsujikawa, Phys. Rev. D 
71, 063004 (2005)]. We show that for a > — 1 there is no Big Rip singularity although w < — 1. For 
n = — 1, we provide an analytical model of phantom bouncing universe "disappearing" at t = 0. We 
also determine the potential of the phantom scalar field and phantom tachyon field corresponding 



to the generalized equation of state p — (ap + kp 

I. INTRODUCTION 

In previous papers of this series, we have constructed 
models of universe with a generalized equation of state 



1 + 1/ r. 



p = (ap + kp 



1+1/™ 



(1) 



having a linear component and a polytropic component. 
In Papers I and II, we have assumed a + 1 + kp 1 / 11 > 
corresponding to w = p/pc 2 > — 1. In that case, the 
density decreases as the universe expands. For n > 0, 
the polytropic component dominates the linear compo- 
nent when the density is high: This describes the early 
universe (Paper I). For n < 0, the polytropic compo- 
nent dominates the linear component when the density 
is low: This describes the late universe (Paper II). When 
the polytropic pressure is positive (k > 0), the solutions 
of the Friedmann equations exhibit past or future singu- 
larities (or peculiarities). When the polytropic pressure 
is negative (k < 0), there is no singularity. Further- 
more, the polytropic equation of state implies the exis- 
tence of an upper bound p m ax (in the past) and a lower 
bound p m i n (in the future) for the density. It makes sense 
to identify the maximum density to the Planck density 
Pp = 5.16 10 99 g/m 3 and the minimum density to the 
cosmological density pa = 7.02 10 _24 g/m 3 . These con- 
stant densities imply in turn the existence of two phases 
of exponential inflation, one in the early universe and one 
in the late universe. During the inflation, the universe is 
accelerating. The early inflation is necessary to solve no- 
torious difficulties such as the singularity problem, the 
flatness problem, and the horizon problem [l|, The 
late inflation is necessary to account for the observed ac- 
celerating expansion of our universe Q driven by dark 
energy [J]. In that context, the equation of state ([T]) 
with k < and n < corresponds to the generalized 
Chaplygin gas Q that has been proposed as a model for 
dark energy. From the generalized polytropic equation 
of state ((T|) , we have obtained a model of universe with- 



out singularity that possesses striking "symmetries" be- 
tween the past and the future (aioniotic universe). This 
model, which could have been obtained from a principle 
of "simplicity" without making any observation, turns 
out to be strikingly consistent with what we know of the 
real universe. It is consistent with the standard model 
@j 0] but refines it by removing the primordial singular- 
ity (Big Bang). In this model, the Planck density and 
the cosmological density are interpreted as two funda- 
mental bounds for the density determined by the Planck 
constant Ti (microphysics) and the cosmological constant 
A (cosmophysics) , respectively. These bounds differ by 
122 orders of magnitudes, a difference that appears to be 
quite natural instead of representing a "problem" Q . 

In this paper, we consider a case that has not been 
treated in our previous papers. This is the case where 
the density increases as the universe expands. Since the 
nature of dark energy is unknown, this situation can- 
not be rejected a priori. It corresponds to an equation 
of state parameter w less than —1 which violates the 
null dominant energy condition. This is referred to as 
a "phantom universe" 0] because when the equation of 
state with w < — 1 is constructed in terms of a scalar 
field, the corresponding kinetic term has the wrong sign 
(negative kinetic energy) . It represents therefore a phan- 
tom (ghost) scalar field (see reviews [il. [Toj). 

Actually, there is a rich recent literature on this sit- 
uation (more than one thousand papers are related to 
phantom dark energy) since observations do not exclude 
the possibility that we live in a phantom universe. In- 
deed, observational data indicate that the equation of 
state parameter w lies in a narrow strip around w = — 1 
possibly being below this value [ll|. The models based 
on phantom dark energy usually predict a future singu- 
larity in which the scale factor, the energy density, and 
the pressure of the universe become infinite in a finite 
time. This would lead to the death of the universe in a 
singularity called "Big Smash" [l2j, "Big Rip" or "Cos- 
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mic Doomsday" [13[ . Contrary to the "Big Crunch" , the 
universe is destroyed not by excessive contraction but 
rather by excessive expansion. In phantom cosmology, 
every gravitationally bound system {e.g. the solar sys- 
tem, the Milky Way, the local group, galaxy clusters) is 
dissociated before the singularity Uj] , and the black 
holes gradually lose their mass and finally vanish [IH, [l6| . 
This scenario allows the explicit calculation of the rest of 
the lifetime of our universe. Actually, as we approach the 
singularity, the energy scale may grow up to the Planck 
one, giving rise to a second quantum gravity era. Eventu- 
ally, quantum effects may moderate or even prevent the 
singularity [It}. Other aspects of phantom cosmology 
have been studied in fl8l ]. 

There are many interesting recent works on the study 
of singularities. In particular, Nojiri et al. fl9j consid- 
ered an equation of state of the form p = —p — f(p) and 
obtained a classification of finite-time future singularities 
(see complements in poj). They are of four types: 

• Type (Big Bang or Big Crunch): For t — > t s , a — > 0, 
p — > +oo, and \p\ — > +00. 

• Type I (Big Rip): For t — > t s , a — > +00, p 
and \p\ +00. 

• Type II (sudden singularity): For t — > t s , a - 
p — > p s , and \p\ — > +00. 

• Type III (Big Freeze): For t — > t s , a — >• a s , p — >• 
and \p\ — > +00. 

• Type IV (generalized sudden singularity) : For t 



+00, 



DO, 



ts 



p s , \p\ — > p s , and higher derivatives of H 



a —> a 
divergi 

In this classification, t s , a s , p s , and p s are all finite 
constants (a s 7^ 0). Type is the standard Big Bang 
or Big Crunch singularity arising in the original Fried- 
mann models @ . Type I is the Big Rip singularity which 
emerges from the phantom equation of state p = ape 2 
with constant a < — 1 [H, [l3[ , and from the equation of 
state (HJ) with a = -1, k < and n < — 2 [2lj. Type 
II corresponds to the sudden future singularity found by 
Barrow [22J at which a and p are finite but p diverges. 
Typ e III, arising in the equation of state (UJ with n > 
[2lL [23( 1 differs from the sudden future singularity in the 
sense that p diverges. Type IV appears in the model 
described in [l9[ . 

It is important to stress that the phantom models with 
W < — 1 do not necessarily lead to future singularities. 
For example, the equation of state ([T]) with a = — 1 and 
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— 2 < n < does not present future singularity 
However, since the scale factor and the density increase 
indefinitely, this has been called "Little Rip" [241 ] . 

On the other hand, the models with w > — 1 may lead 
to past or future singularities. For example, the new 
form of primordial singularity (for n > and k > 0) 
described in Sees. IV C, VI and in Appendix A of Paper I 
corresponds to a past singularity of type III: The universe 



1 We shall not consider this type of singularities in this paper. 



starts at t = with a finite scale factor and an infinite 
density. On the other hand, the future singularity (for 

— 1 < n < and k > 0) described in Sec. IV C and in 
Appendix B of Paper II corresponds to a singularity of 
type II: At a finite time t s , the universe reaches a point 
at which the scale factor is finite, the density vanishes 
and the pressure is infinite. 

In paper II, we have also introduced a notion of "pecu- 
liarity" . This is when the density vanishes p — while the 
scale factor has a finite value a s (when a s = we shall call 
it generalized peculiarity). In that case, the universe is 
empty (in other works, it "disappears"). Although there 
is no singularity, this situation is very peculiar. However, 
since the nature of dark energy is unknown, all possibil- 
ities should be contemplated. 

In this paper, we provide an exhaustive study of the 
equation of state (fTJ in the case w < — 1 (requiring k < 
0) for arbitrary — 1 < a < 1 and n. This is a natural 
continuation of our previous works which assumed w > 

— 1 (Papers I and II). Our paper also completes previous 
studies of the case w > — 1 that considered a = — 1 
[2l| or a = [23|. An interesting result of our study 
is that the equation of state (fTJ) with a > — 1 does not 
present a Big Rip singularity although w < — 1, contrary 
to the linear equation of state p = ape 2 with a < — 1 
fl3| or the equation of state (fTJ) with a = — 1 and n < 
—2 2l|. Another interesting result of our study is the 
construction of a bouncing phantom universe for — 2 < 
n < 0. For — 1 < n < 0, this bouncing universe presents 
a past singularity of type II: At t = 0, the the pressure is 
infinite while the scale factor has a finite value and the 
density vanishes. For a > — 1 and n = —1, corresponding 
to a constant negative pressure, the bouncing phantom 
universe admits a simple analytical expression. 

The paper is organized as follows. In Sec. HH we recall 
the basic equations of cosmology. In Sees. IIIII and IIV1 
we study the generalized equation of state (fTJ) for any 
value of the parameters — l<a<l,fc<0 and n, 
assuming w < — 1 (phantom cosmology). In Sec. El we 
determine the potential of the phantom scalar field and 
the potential of the phantom tachyon field corresponding 
to the generalized equation of state (fTJ) . In Appendix 
we treat the case a = — 1. In Appendix [B] we summarize 
all the results obtained in our series of papers and analyze 
the different singularities in terms of the classification of 

El. 



II. BASIC EQUATIONS OF COSMOLOGY 

We assume that the universe is isotropic and homo- 
geneous at large scales and contains a uniform perfect 
fluid of energy density e(t) = p{t)c 2 and pressure p(t). 
We also assume that the universe is flat in agreement 
with observations of the cosmic microwave background 
(CMB) 0. Finally, in this paper, we ignore the cos- 
mological constant (A = 0). It that case, the Einstein 
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equations reduce to 



— + 3- ( + — 

dt a V c 2 



= 0, 



H- 



AttG 



8ttG 



-P, 



(2) 
(3) 

(4) 



where a(t) is the scale factor ("radius" of the universe) 
and H = a /a is the Hubble parameter. These are the 
well-known Friedmann equations describing a non-static 
universe @. The first equation can be viewed as an 
"equation of continuity" . For a given barotropic equa- 
tion of state p = p{p), it determines the relation between 
the density and the scale factor. Then, the temporal evo- 
lution of the scale factor is given by Eq. (Q| . Introducing 
the equation of state parameter w = p/ pc 2 , we see from 
Eq. ([3]) that the universe is decelerating if w > —1/3 
(strong energy condition) and accelerating if w < —1/3. 
On the other hand, according to Eq. @, the density 
decreases with the scale factor if w > — 1 (null dominant 
energy condition) and increases with the scale factor if 
w < — 1. In this last case, we are dealing with a "phan- 
tom universe" . 

We will also need the thermodynamical equation 



1 



dp 



(pc 2 +p), 



(5) 



which can be derived from the first principle of thermo- 
dynamics For a given barotropic equation of state 
p = p(p), this equation can be integrated to obtain the 
relation T = T(p) between the temperature and the den- 
sity. It can be shown Q that the Friedmann equations 
conserve the entropy of the universe 



S= a Y {p + pc% 



(6) 



If we impose that the entropy is positive, we conclude 
from Eq. ([6]) that the temperature is positive when 
w > — 1 while it is negative when w < — 1. The fact 
that a phantom universe has a negative temperature was 
mentioned in Negative temperatures arise in other 
domains of physics such as 2D turbulence [25| . 

The simplest model of phantom universe corresponds 
to the linear equation of state p = ape 2 with a < — 1 
The continuity equation ([2]) can be integrated into 



p oc a 



3\l+a\ 



(7) 



Substituting Eq. (J7J in Eq. |QJ and solving the resulting 
equation for a(t), we find that the scale factor, the Hubble 
parameter and the density increase in time as 



a oc (t s — t) 



-2/(3|l+o| 



(8) 



H 



P 



3 1 + a 



(ts-t)- 1 , 



6ttG(1 + a-) : 



(9) 



(10) 



They all diverge at a finite time t = t s . This is the "Big 
Rip" singularity [lH , which is a singularity of type I [l|| . 
We also find from Eq. ([5} that the temperature behaves 
as 



Tex 



a/(a+l) 



oc -a 3|Q| OC -(t s -t) 



-2a/(l+a) 



(11) 



The temperature becomes more and more negative as the 
universe expands, and it diverges when t—>t s . 



III. GENERALIZED EQUATION OF STATE 
WITH w< -1 

We consider a generalized equation of state of the form 

p= {ap + kp 1+1 ' n )c 2 . (12) 

This is the sum of a standard linear equation of state 
p = ape 2 and a polytropic equation of state p = kp^c 2 , 
where k is the polytropic constant and 7 = 1 + 1/n is 
the polytropic index. Concerning the linear equation of 
state, we assume — 1 < a < 1 (the case a = — 1 is treated 
specifically in Appendix [5J ■ Concerning the polytropic 
equation of state, we remain very general, so that k and n 
can take arbitrary values. In papers I and II, we assumed 
that a + 1 + kp x l n > 0, so that the density decreases 
with the scale factor (w > — 1). In the present paper, we 
assume that a + 1 + fcp 1 /™ < (a necessary condition is 
k < 0) so that the density increases with the scale factor 
(w < —1). This corresponds to a "phantom universe". 

A. The density 



For the equation of state (|T2|) , the Friedmann equation 
@ becomes 



^ + ^p(l + a + kp^ n ) = 0. 
dt a 



(13) 



Assuming a + 1 + kp 1 ! 11 < 0, this equation can be inte- 
grated into 



P 



[1 - (o/o.) 3 C 1 +»)/»] 



(14) 



where = [(a+ l)/|fc|]™ and a* is a constant of integra- 
tion. 

For n > 0, the density is defined only when a < a*. 
When a — > 0, p — > p* and p — > —p*c 2 . When a — > a*, 



P_ 

P* 



1 



.3(1 + a). 



(1 - a/a«) n 



(15) 



a/a, 

FIG. 1. Density as a function of the scale factor for n > 
and n < (specifically n = 1 and n = —1/2). We have taken 
a = 0. 



a/a, 

FIG. 2. Temperature as a function of the scale factor for 
n > and n < (specifically n = 1 and n = —1/2). We have 
taken a = 0. 



and p — > — oo. 

For n < 0, the density is defined only when a > a*. 
When a — >• a*, 



/>* 



3(1 + a) 



M 



(a/a* - l) |n| -> 0. 



(16) 



In the same limit, p — > — oo for n > — 1, p tends to a 
finite value for n = — 1, and p — > for n < — 1. On the 
other hand, when a — > +oo, p — > p* and p — > — p*c 2 . 

Some curves giving the evolution of the density p as 
a function of the scale factor a are plotted in Fig. [1] for 
n > and n < 0. 



B. The temperature 



For the equation of state (|12[) . the thermodynamical 
equation ([5]) can be integrated into 



T = -T* 



a/n _ 



(a+n+l)/(a+l) 



(P/P*) 



q/(q+1) 



(17) 

where J 1 * > is a constant of integration. Combined 
with Eq. (|14p . we obtain 



r 



-t. 



(a/a, 



\3{a+n+l)/n 



[1 - (a/a,)3(i+«)A 



1 n+l 



(18) 



We have to consider different cases. 

We first assume n > 0. When a — > 0, T — > 0; when 
a — >• a* , T — > — oo. 

We now assume n < 0. When a — > a*, T — > for 
n < —1 and T — > — oo for n > — 1. When a — > +oo, 
T^0forn + Q + l>0 and T -> -oo for n + a + 1 < 0. 

The extremum of temperature (when it exists) is lo- 
cated at 



Pe 
P* 



(l + a)(n + l) 



(19) 



a P I a + n + 1 

a* 



an 



n/[3(l+a)] 



(20) 



Te 



n + a + 1 



i+a+l _ 
1+a 



(a + l)(n + l) 



n+l 



(21) 



Some curves giving the evolution of the temperature T 
as a function of the scale factor a are plotted in Fig. [2] 
for n > and n < 0. 

Finally, the entropy ([5]) is given by 



(22) 



and we explicitly check that it is a positive constant. 



C. The parameter w(t) 

We can rewrite the equation of state (fi"2j) as p = 
w(t)pc 2 with 



w(t) = a - (a + 1) 



l/n 



(23) 



For n > 0, w — > — 1 when a — > and w — > — oo when 
a — > a*. 

For n < 0, it) — > — oo when a — > a* and it) — > —1 when 
a — > +oo. 

Some curves giving the evolution of w as a function of 
the scale factor a arc plotted in Figs. [3] and 0] for n > 
and n < 0. 
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-(a+n+l)/n 1 \\ 









0.5 



a/a,,, 



FIG. 3. The parameters w, q and c? a /c? as a function of the 
scale factor a for n > (specifically n = 1). We have taken 
a = 0. 
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Vjj< /C) -(a+n+l)/n 
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r. , ■ 
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a/a. 



FIG. 4. The parameters to, g and c? 8 /c? as a function of the 
scale factor a for n < (specifically ?i = —1/2). We have 
taken a — 0. 



D. The velocity of sound 

For the equation of state (fT2"j) , the velocity of sound is 
given by 



P'(p) 



a - (a + 1) 



n + 1 



1/n' 



(24) 



For n < 0, the velocity of sound vanishes at the point 
(fT9")) - (|21[) where the temperature is extremum. At that 
point, the pressure is maximum with value 



Pe 

p*c 2 



71+1 



(l + a)(n+l) 



(25) 



The case c 2 < corresponds to an imaginary velocity of 
sound. We also define 



Ps 

p* 



(1 — a)n 



(l + a)(n + l) 



(26) 



a + 2n + 1 



a) 



n/[3(l+a)] 



(27) 



corresponding to a possible point where the velocity of 
sound is equal to the speed of light. Different cases have 
to be considered. 

We first assume n > 0. When a — > 0, (c s /c) 2 — > — (a + 
n + l)/n; when a — > a*, (c s /c) 2 — > — oo. The velocity of 
sound is always imaginary. 

We now assume n < 0. When a —» a*, (c s /c) 2 — ^ 
+co for n > —1 and (c s /c) 2 — > — oo for n < — 1; when 
a —> +oo, (c s /c) 2 — > —(a + n + For n > — 1 and 

a + n + 1 > 0, c 2 in always positive. For n > — 1 and 
a + n + 1 < 0, it is positive for a < a e and negative for 
a > a e . For n < —1 and a + n + 1 < 0, c 2 in always 
negative. For n < —1 and a + n + 1 > 0, it is negative 
for a < a e and positive for a > a e . For ri > — 1 and 
a + In + 1 > 0, the velocity of sound is always larger 
than the speed of light. For n > — 1 and a + 2n + 1 < 0, 
the velocity of sound is larger than the speed of light 
for a < a s and smaller for a > a s . For n < — 1 and 
a + 2n + 1 < 0, velocity of sound is always smaller than 
the speed of light. For n < — 1 and a + 2n + 1 > 0, the 
velocity of sound is smaller than the speed of light for 
a < a s and larger for a > a s . 

Some curves giving the evolution of (c s /c) 2 as a func- 
tion of the scale factor a are plotted in Figs. [3] and 0] for 
n > and n < 0. 



IV. EVOLUTION OF THE SCALE FACTOR 

A. The deceleration parameter 

The deceleration parameter is defined by Eqs. (1-77) 
and (F78). A phantom universe is always accelerating 
since q < — 1 < 0. For the equation of state (fT2"]l . using 
Eq. we get 



1 + 3a 3 . , / p 

?(* ) = — 5 n "+ 1 - 

2 2 \p* 



l/r 



(28) 



For n > 0, q — > —1 when a — > and g — > — oo when 
a — > a*. 

For n < 0, g — > — oo when a — > a* and q — > — 1 when 
a — > +oo. 

Some curves giving the evolution of q as a function of 
the scale factor a are plotted in Figs. [3] and 0] for n > 
and n < 0. 



B. The differential equation 

The temporal evolution of the scale factor a(t) is de- 
termined by the Fricdmann equation (U). Introducing 
the normalized radius R = a/a*, the density (|14p can be 
written 

P* 



P = 



[1 - i? 3 (!+ Q )A 



(29) 
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Substituting this expression in Eq. we obtain the 

differential equation 

eKR 



R 



(30) 



[1 - j R3(l+«)/n]n/2' 

where K = (SnGp^/S) 1 ^ 2 and e = ±1. In general, we 
shall select the sign e = +1 corresponding to an expand- 
ing universe (R > 0), except in the case of a bouncing 
universe where both signs of e must be considered. The 
solution can be written as 

V2 dR 



eKt 



1 — i? 3 ( 1+Q )/ n 



R ' 



or, after a change of variables x — i? 3 ( 1+Q )/", as 



3(a + i; 



eKt 



(1-x) 



re/2 



dx 



(31) 



(32) 



The integral can be expressed in terms of hypergeomet- 
ric functions. Some simple analytical expressions can be 
obtained for specific values of n. Actually, we can have a 
good idea of the behavior of the solution of Eq. f|30|> by 
considering asymptotic limits (see below). The complete 
solution is represented in the figures by solving Eq. (|30p 
numerically. 



C. The case n > 

The universe starts from t — > — oo with a vanishing 
radius R = 0, a finite density p = p*, and a finite pressure 
p = —p*c 2 . When t — > — oo, 

R ~ Ae Kt . (33) 

This corresponds to an exponential expansion (early in- 
flation) due to the fact that the density is approximately 
constant. Then, the universe undergoes a finite time sin- 
gularity at a time When t — > £*, the radius tends to 
its maximum value R = 1 while the density tends to +oo 
and the pressure to — oo. This is a future singularity of 
type III. Close to the singularity, we have 

n 2/(2+re) 



R 



p* 



3(a + l) 



-(1 + a) 

2 n 



i/2 



K(t* - 1) 



'■K(U-t) 



-2ra/(2+n) 



. (34) 



(35) 



The evolution of the scale factor is represented in Fig. 
[SJ Some simple analytical results can be obtained in par- 
ticular cases. 

For n = 1, using the identity 

■ dx 



Vl-x — = 2V1 - x 

x 



In 



i - yi^ 

l + v 7 ! - 



and for n = 2, using the identity 

dx 



(1-x) 



-x + lax, 



(36) 



(37) 



we can obtain t(R) from Eq. (f3"2")l . 



l 

0.8 - 
[*!0.6 
0.4 
0.2 




n>0 



-2 
Kt 



FIG. 5. Evolution of the radius R as a function of time for 
n > (specifically n = 1). We have taken a = 0. We 
have chosen the constant of integration such that R = 1 at 
t = U = l/K. 



D. The case n < 

The early evolution of the universe depends on the 
value of n. Different cases must be considered. 

(i) For n < —2, the universe starts from ( 4 -co 
with a finite radius R = 1, a vanishing density p = 0, 
and a vanishing pressure p = (past peculiarity) . When 



t 



-oo, 



R-l 



3(a + l) 



i\/2 



2/(M-2) 



(-Kt) 



P* 



1 



1 



3a + 1 \n\ - 2 (-ift) 



2|n|/(M-2) 



(38) 



(39) 



The density tends to zero algebraically rapidly. 

(ii) For n = —2, the universe starts from t — > — oo 
with a finite radius R = 1, a vanishing density p = 0, 
and a vanishing pressure p — (past peculiarity) . When 
i — » — oo, 



R - 1 ~ 



-(l + a) 2 ^^ 3 ^)^. 



(40) 



(41) 



The density tends to zero exponentially rapidly. 

(hi) For — 2 < n < 0, the universe starts at t = with 
a finite radius i? = 1 and a vanishing density p = (past 
peculiarity). When t — > 0, 



i?- 1 



2-N 



3(a + l) 



|re|/2 1 2/(2-|n|) 

A'i I , (42) 



3, .2-|n| T , ' 
_(l + a —Uii-t 

2 n 



2|n|/(2-|n|) 



(43) 
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At t = 0, the pressure vanishes for — 2 < n < — 1, is finite 
for n = — 1 and tends to — oo for n > — 1. In this last 
case, there is a past singularity of type II. Actually, we 
can extend the solution to t < (except, maybe, in the 
case n > — 1 where the pressure diverges). This describes 
a phase of contraction of the universe, corresponding to 
the solution of Eq. (PU|) with e = — 1. This leads to a 
model of bouncing phantom universe that collapses for 
t < (with decreasing density), disappears at t = (the 
density vanishes), and expands for t > (with increasing 
density). 

On the other hand, for t — > +oo, the density tends to a 
constant p* , implying an exponential growth of the scale 
factor as 

R ~ A!e Kt . (44) 

This corresponds to a phase of late inflation. The pres- 
sure p — > —p*c 2 . 




Eq. (|46)) has the asymptotic forms (|42|) and (|44|) with 
A! = 2~ 2 /[ 3 ( 1+Q )l. For a = 0, this model has a constant 
negative pressure p = — \k\c 2 . It belongs therefore to 
the same "class" as the ACDM model (Paper I) and the 
anti-ACDM model (Paper II) that also have a constant 
pressure. These three models admit simple analytical 
expressions. 




Kt 



FIG. 7. Analytical model of a bouncing phantom universe 
corresponding to n = — 1. We have taken a = 0. The universe 
starts from t = — oo with a maximum density p*. It first 
experiences a phase of contraction during which its density 
decreases. At t = 0, it reaches its minimum radius R = 1 
and its density vanishes (the universe "disappears"). For t > 
0, the universe expands and its density increases up to the 
maximum value p, . 



Kt 



For n = — 2, using the identity 



FIG. 6. Evolution of the radius R as a function of time for 
n < —2, n — —2 and —2 < n < (specifically n = —3, 
n = —2 and n = —1/2). We have taken a = 0. For n < —2, 
we have chosen the constant of integration such that R — 2 
at t = 0. For —2 < n < 0, the constant of integration has 
been chosen such that R = 1 at t = 0. In that case, we have 
a bouncing universe (see text for details). 



The evolution of the scale factor is represented in Fig. 
|5] for n < 0. Some simple analytical results can be ob- 
tained in particular cases. 

For n = — 1, using the identity 



1 



dx 



= in 



we obtain 



2/[3(l+a)] 



i? = cosh 



— = tanh 2 
P* 



l + VT 



■(1 + a)Kt 



-(l + a)Kt 



(46) 



(47) 



This provides an analytical solution of a bouncing phan- 
tom universe (see Fig. [7]). We can explicitly check that 



1 dx 

m 



we obtain 



1 — x x 



r= i + e i(i+«)** 



P* 



1-x ' 



2/[3(1+q)] 



1 



\{l+a)Kt 



(48) 



(49) 



(50) 



This provides an analytical solution of a phantom uni- 
verse exhibiting a past peculiarity and a late inflation. 
Equation (|49| has the asymptotic forms (|40)) and (|44)) 



(45) with A' = 1 and A = 2/[3(l + a)]. 



For n = —1/2, we have the identity 



1 



(l-x) 1 /* x 



— = 2 arctan 



-In 



(1-x) 1 / 4 

1 - (1-x) 1 / 4 
l + (l-x)V 4 



(51) 



which determines t(R) using Eq. (|32|) . For a = 0, this 
solution corresponds to the phantom Chaplygin gas. 

Remark: We note that the solution with n < — 2 looks 
similar to the Eddington-Lcmaitrc model (see Fig. 1 in 
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Paper I) since the universe is "static" in the past with a 
finite radius R = 1 and expands exponentially rapidly in 
the future. However, in the Eddington-Lcmaitrc model, 
the density decreases with time while, in the present 
(phantom) model, it increases with time. In addition, 
in the Eddington-Lemaitre model, the density tends to 
a finite value when t — > — oo while in the present model, 
it tends to zero. Therefore, these models are physically 
very different. 



V. SCALAR FIELD MODELS 

In this section, we introduce a representation of the 
phantom universe in terms of scalar field models. We de- 
termine the potential of the scalar field corresponding to 
the equation of state (TT2"]) using the general methodology 
exposed in Q. We consider a normal scalar field and a 
tachyon field. Although totally equivalent to fluid equa- 
tions, this scalar field representation may be useful in 
order to make the link with more fundamental theories, 
like those arising in particle physics and string theory 
0. 



A. Phantom scalar field 

A fluid with an equation of state parameter satisfying 
w > — 1 can be described in terms of an ordinary scalar 
field minimally coupled to gravity called a quintessence 
field [26|. A fluid with an equation of state parameter 
satisfying w < — I can be described in terms of a phantom 
scalar field. The phantom field can be obtained from the 
quintessence field by making the transformation <fi —> i<p. 
As a result, the phantom scalar field evolves according to 
the equation 



3^-^ = 0, 



(52) 



where V(4>) is the potential of the scalar field. The den- 
sity and the pressure are given by 



pc 2 = -l<j> 2 + V(<f>), 



P = -U 2 -Vttl ( 53 ) 



We note that the sign of the kinetic term is reversed with 
respect to the quintessence scalar field. This implies that 
the phantom scalar field tends to run up, not down, the 
potential towards larger energies. 
From Eq. (|53[) . we get 



cj> 2 = |1 + w\pc 2 



(54) 



where we have written p = wpc 2 . Using tj> = (d(f>/ 'da)Ha, 
and the Fricdmann equation Q valid for a flat universe, 
we obtain 



da 



3c 2 
SttG 



2 \ 1/2 



(55) 



For the equation of state ([12]), using Eqs. ([14]) and ([23]) . 
and setting R = a/a*, we can rewrite Eq. (|55p in the 
form 



3c 2 \ 1/2 V^TT i?3(l+«)/2« 



dR \8ttG J R v 7 ! - i?3(l+a)/n 

With the change of variables 

1/2 



(56) 



^(a+D/an J, = ^±±0, (57) 

\ 3c 2 J 2n 



we find that 



1> 



dx 



VI - x< 



= Arcsin(x) 



On the other hand, according to Eq. ([53]) . we have 

V = l -(l-w)pc 2 . 

For the equation of state lfT2"]) . using Eqs. (fT4"]l and 
we obtain 

1 2 2-(l-a)x 2 

2 P * C (I- X 2 )^ 1 • 



(58) 



(59) 



(60) 



Since x = sin ip, the scalar field potential is explicitly 
given by 



I 2 (I — a) cos 2 ip + a + 1 



cog 2(n+l) ^ 



(61) 



and ij3(o+i)/2n = sin ^,^ In thesc mo dcls, < -0 < tt/2. 

The case a = — 1 and k < (sec Appendix |A]) must be 
treated specifically. Repeating the preceding procedure, 
we find that the potential of the scalar field is 



v{4>) = P *c 2 



\n\c 2 
2^G 



2 \ ™ 



2^,2 



n c 



2 / A2- 



lni?= — -(f) 2 . 

nc 2 



(63) 



In these models (j> > 0. 

Finally, for a linear equation of state p = ape 2 with 
a < —I, writing the relation between the density and 
the scale factor as p/p* = (a/a*) 3 ' 1+a l, we obtain Q: 

= I^ c 2 (1 _ a )e^V^\m) 1/2 ^ (64) 



8itGJ 



\nR, 



(65) 



where R = a/a* and <j> < 0. Since R cx (t s _f)-2/(3|i+a|) ; 
the scalar field evolves in time as <j> = — (c 2 /6irG\l + 
a]) 1 / 2 ln(t s -t). 
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B. Phantom tachyon field 

Performing the transformation <f> — > i<f> in the equations 
of an ordinary tachyon field p7j , we find that a phantom 
tachyon field evolves according to the equation 



The density and the pressure are given by 



pc 



P 



! l + (f> 2 

From these equations, we obtain 

4> 2 = + H 



v{4>)di + 4> 2 - 



(66) 



(67) 



(68) 



where we have written p = wpc 2 . Using 4> = (d<f>/da)Ha, 
and the Friedmann equation (U]), we get 



da 



3c 2 
SttG 



2 \ 1/2 



' pc 2 a 



(69) 



For the equation of state (IT21) . using Eqs. (fT¥|) and 
we can rewrite Eq. (|69[) in the form 



d<j> 



3c 2 ^ 1/2 



+ \ p3(l+q)/2n 



X 1 - J R 3 ( 1+Q )/" 

With the change of variables 



(n-l)/2 



(70) 



0= / (l-z 2 )^" 1 )/ 2 ^. 



we find that 



On the other hand, from Eq. (|57)) . we have 
U 2 = - W p 2 c 4 . 



(71) 
(72) 

(73) 



For the equation of state (H2]), using Eqs. (HU) and (|2"3"|) . 
we obtain 



2„4 



oke + 1 



(74) 



Therefore, the scalar field potential V (i/j) is given in para- 
metric form by Eqs. (|72|) and (|74|) . Let us consider par- 
ticular cases. 

(i) For n = 1, we find that x = i\). Therefore, we obtain 



U 2 = p 2 c 4 °^ + 1 



(i -V 2 ) 3 ' 



(75) 



and i?3(«+i)/2 _ ^ with < V < 1. 

(ii) For n = — 1, we find that x = tanh-0. Therefore, 
we obtain 



2 4 a tanh -0 + 1 



cosh ^ 



(76) 



and i?3(a+i)/2 _ i/tanhip with > 0. 

(iii) For n = — 2, we find that a; = ip/yl + iJj 2 . There- 
fore, we obtain 



2 2 4 (a + 1)^ + 1 

P * (1+^2)4 > 



(77) 



and i?3(a+i)/4 = ^1 + ^2/^ with ^ > o. 

(iv) For n = —1/2 and a = (phantom Chaplygin 
gas), we find that V{<j>) — p*c 2 is constant. 

The case a = — 1 and fc > (see Appendix IX)) must be 
treated specifically. Repeating the preceding procedure, 
we find that the potential of the scalar field is 



x 1 + T 



2irGp«(n + l) 2 
H 

2nGp4n + l) 2 



2n/(ra+l) 



l/(n+l) 



n+1 



■]nR = sgn(n) 



2irGp*(n+ if 



l/(n+X) 



A2/("+l) 



(78) 



(79) 



In these models </> > 0. 

Finally, for a linear equation of state p = ape 2 with 
a < — 1, writing the relation between the density and 
the scale factor as p/p* = (a/a*) 3 l 1+Q l, we obtain j4[: 



2 1 



y(0) = 



3c 2 



c 2 1 
1 + a\ 6ttG <j> 2 ' 



(80) 



^ =jR -3|i+ Q |/2 ; (81) 



where R = a/a* and 4> < 0. Since p = p*i? 3 l 1 + Q l = 
l/[67rG(l + a) 2 ](t s — 1)~ 2 , the scalar field evolves in time 
as 4> = -«J\l + a\ (t s -t). 



VI. CONCLUSION 

In this paper, we have performed an exhaustive study 
of the generalized equation of state (fT"2|) in the case where 
the pressure increases with the scale factor. This corre- 
sponds to the so-called phantom cosmology @. 

The case a = — 1 was previously treated in [2l|. For 
n < — 2, the universe experiences a future singularity of 
type I (Big Rip): The scale factor and the density diverge 
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at a finite time. For —2 < n < 0, the scale factor and the 
density diverge in infinite time (Little Rip). For n > 0, 
the universe experiences a future singularity of type III: 
The density diverges at a finite time while the scale factor 
tends to a constant. 

We have found that when a > — 1, the universe does 
not experience a Big Rip singularity. For n < 0, there 
is a phase of late inflation and, for n > 0, the universe 
experiences a future singularity of type III. The past evo- 
lution of the universe is interesting. For n > 0, there is 
a phase of early inflation. For n < 0, the universe ex- 
hibits a past peculiarity since the density vanishes while 
the scale factor tends to a finite value. For n < — 2, the 
evolution of the scale factor is similar to the Eddington- 
Lcmaitre model (the universe is static in the infinite past 
and grows exponentially rapidly in the future) but the 
evolution of the density is very different (it starts from 
zero in the infinite past and increases as the universe ex- 
pands). For —2 < n < 0, we have obtained a model of 
bouncing universe which also possesses peculiar features 
(the density decreases in the past, vanishes at t = 0, and 
increases in the future). For — 1 < n < 0, this bouncing 
universe presents a past singularity of type II since the 
pressure at t = is infinite while the scale factor is finite 
and the density vanishes. For n = — 1, corresponding 
to a constant negative pressure, the bouncing phantom 
universe admits a simple analytical expression. 

Of course, most of these models are academic, and 
do not correspond to the true evolution of our universe. 
However, we believe that it is important to study the 
equation of state (fT2|) in full generality. On the other 
hand, there are indications [ll| that the equation of state 
parameter w of our universe may become less than — 1 in 
the close future (or even at present), so the late evolution 
of the phantom models described in this paper may be 
physically relevant. 

A drawback of the simple form of phantom cosmol- 
ogy considered in this paper is that it does not connect 
smoothly to the matter era (which has w — 0). There- 
fore, we cannot realistically extend the phantom models 
to the past and obtain unified models of dust matter 
and phantom dark energy (w < —1), contrary to the 
unified models of dust matter and quintessence dark en- 
ergy (w > —1) based on the generalized Chaplygin gas 
considered in Paper II. A unification of dust matter and 
phantom dark energy can be achieved in mo re g eneral 
models allowing to cross the phantom divide [28|. This 
generalization assumes an interaction between dark mat- 
ter and dark energy. These models are very interesting 
because they may provide a solution to the "cosmic coin- 
cidence problem" (the fact that the ratio of dark matter 
and dark energy is of order one). 

The phantom cosmology is also interesting for its con- 
nection to Hoyle's version of the steady state theory [2!| , 
for its connection to wormholes [30l ]. and for its very 
strange thermodynamics allowing for the existence of 
negative temperatures [l6| like in 2D turbulence 25 1. 

However, we may recall that there is no firm evidence 



that we live in a phantom universe. The model of Paper 
II, corresponding to the standard ACDM model with the 
primordial singularity removed, may correctly describe 
the whole evolution of our universe. Therefore, a more 
precise determination of the equation of state parameter 
w will help discriminate between these different models. 



Appendix A: Equation of state p = (— p + kp~')c 2 with 

k < 

In this Appendix, we specifically study the equation of 
state (|T2|) with a = — 1 and k < 0, namely 



p = (-p-|fe|/>V 



(Al) 



Since w < — 1 , this equation of state describes a phantom 
universe. This equation of state was introduced by No- 
jiri & Odintsov [TtJ and studied by Stefancic [2l|. Nojiri 
et al. pjl used it to illustrate their classification of fu- 
ture finite time singularities. For the completeness of our 
study, we shall re-derive their results in a more compact 
form (with our notations) and give a few complements. 
We follow the same presentation as in Papers I and II. 



1. The case n > 

The equation of continuity ^ can be integrated into 

P* 



ln(a*/a) n 



(A2) 



where = (n/3|fc|)" and a* is a constant of integration. 
The density is defined for a < a*. When a — > 0, p — > 
and p — > 0; when a — ¥ a*, p — > +oo and p — > — oo. 



H 
H 



1 


1 1 1 1 1 1 ' 


n>0 






k<0 




p'p, / 


a = -l 




T7T, \ | 




1 


,,,,,1.1,, 






0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



FIG. 8. Evolution of the density and temperature as a func- 
tion of the scale factor. We have taken n = 1. 



The thcrmodynamical equation ([5]) can be integrated 
into 



T = -T* ( — 

,P 



, (n+l)/n 



(A3) 
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where T» > is a constant of integration. Combined 
with Eq. (| A2|) . we obtain 



T = — 



71 



ln(a*/a) Tl+1 \a 



(A4) 



When a — > a*, T — > 0; when a — > a*, T — > — oo. The 
evolution of the density and temperature as a function of 
the scale factor is represented in Fig. [8j 

The equation of state can be written as p = wpc 2 with 



w 



3 \P* 



l/r, 



(A5) 



When a — > 0, w — s- — 1; when a — > a*, in — > — oo. 

The deceleration parameter is given by Eqs. (1-77) and 
(1-78). Together with Eq. (|A"5j) . we obtain 



2 \p 



l/n 



(A6) 



When a — > 0, g — > —1; when a — > a*, q 
The velocity of sound is given by 



-oo. 



-§ = -! 



n + 1 /> 



(A7) 



When a — s- 0, (c s /c) 2 — » — 1; when a — s- a*, (c s /c) 2 
— oo. The velocity of sound is always imaginary. 



pW 

p* 



2 + n 



K{U - 1) 



-2n/(2+n) 



(A10) 



The universe starts from t — — oo with a vanishing radius 
R = 0, a vanishing density p = and a vanishing pressure 
p = 0. This corresponds to a generalized past peculiarity. 
It also undergoes a future singularity of type III: At t = 
t#, the density tends to +oo and the pressure tends to 
— oo while the radius reaches its maximum value R = 1. 




FIG. 10. Evolution of the radius with time for n > (specif- 
ically n — 1). We have arbitrarily taken t* = 1. 



The evolution of the scale factor is represented in Fig. 




FIG. 9. Evolution of w, q, and (c s /c) 2 as a function of the 
scale factor a. We have taken n = 1. 



The evolution of w, and (c s /c) 2 as a function of the 
scale factor a is represented in Fig. [SJ 

Setting R = a /a* , the Friedmann equation (j4|) can be 
written 



R 



KR 



(A8) 



(-lni?)™/ 2 ' 
where if = (SttGp^/S) 1 / 2 . Its solution is 

R(t) = e-i^ K ^r /i2+ "\ (A9) 



2. The case n < 

The equation of continuity ^ can be integrated into 

p= 1 /; w (ah) 

m(a/a») n 

where p* = (|n|/3|fc|) n and a* is a constant of integration. 
The density is defined for a > a*. When a — >■ a*, p — > 0. 
In the same limit, p — > for n < — 1, p tends to a finite 
value for n = — 1, and p — > — oo for n > — 1. When 
a — > +oo, p — > +oo, and p — > — oo. 

The thcrmodynamical equation ([5]) can be integrated 
into 



T = -T* ( — 

vP 



(n+l)/n 
\ ^(p./p) 1 /' 



(A12) 



where > is a constant of integration. Combined 
with Eq. (fATTl) . we obtain 



T 



ln(a/a*) ,l+1 \a* 



(A13) 



When a — > a*, T — > for n < —1 and T — > — oo for 
n > —1. When a — > +oo, T — > — oo. For n < —1, the 
temperature reaches its maximum at 



Pe 
P* 



n+l 



— = P ("+ 1 )/ 3 



(A14) 
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71+ 1 



n+l 



„n+l 



(A15) 



The evolution of the density and temperature as a func- 
tion of the scale factor is represented in Fig. [TTJ 
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FIG. 11. Evolution of the density and temperature as a func- 
tion of the scale factor. We have taken n = —1/2. 



We have to distinguish several cases. We first assume 
n < — 1. When a — > a*, (c s /c) 2 — > — oo; when a — > +oo, 
(c s /c) 2 —> —1. The velocity of sound is always imaginary. 
We now assume n > — 1. When a — > a*, (c s /c) 2 — > +oo; 
when a — > +oo, (c s /c) 2 — > — 1. The velocity of sound 
vanishes at the point (|A14j) at which the temperature is 
maximum. At that point, the pressure is maximum with 
value 



Pe 

p*c 2 



1) 



n+l ' 



(A19) 



The velocity of sound is real for a < a e and imaginary 
for a > a e . On the other hand, the velocity of sound is 
equal to the speed of light at 



Ps 

p* 



6 



1 



a* 



a (n+l)/6 



(A20) 



The velocity of sound is larger than the speed of light 
when a < a s and smaller when a > a s . The evolution 
of w, q, and (c s /c) 2 as a function of the scale factor a is 
represented in Fig. [T2] 



The equation of state can be written as p = wpc 2 with 

l/n 



w 



-1 



n I p_ 

3 VP* 



(A16) 



When a — > a*, w — > — oo; when a — > +oo, iu — > —1. 

The deceleration parameter is given by Eqs. (1-77) and 
(1-78). Together with Eq. (TAT6]) . we obtain 



2 \P* 



l/n 



(A17) 



When a — > a*, q — > — oo; when a — > +oo, q — > — 1. 



1 1 1 
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n>-l n<0 
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FIG. 12. Evolution of w, q, and (c s /c) 2 as a function of the 
scale factor a. We have taken n = —1. 



k<0 

n<0 




FIG. 13. Evolution of the radius with time for n < 0. We 
have taken n = —3, n — —2 and n = —1/2. For n < —2, we 
have taken Kt* = 1. For n > —2, the universe is bouncing at 
t — (see text for details). 



Setting R = a/a*, the Friedmann equation ([4|) can be 
written 



1? 



eKR 
(lni?)™/ 2: 



(A21) 



where A' = (SttG/^/S) 1 / 2 and e = ±1. We must distin- 
guish three cases, 
(i) For n < -2, 



R(t) = el^^C*--*) 



(A22) 



The velocity of sound is given by 

,2 



-§ = -! 



n + 1 /_p 



l/n 



(A18) 



p* 



Inl -2 



K(U-t) 



-2|n|/(|n|-2) 



(A23) 



The universe starts from t = — oo with a finite radius R = 
1, a vanishing density, and a vanishing pressure (past 
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peculiarity). It undergoes a future singularity of type I 
(Big Rip): At t = t*, the scale factor, the density and 
the pressure are infinite. The divergence of the density 
is algebraic, 
(ii) For n = -2, 



R(t) 



P* 



(A24) 



(A25) 



The universe starts from t = — oo with a finite radius 
R = 1, a vanishing density, and a vanishing pressure 
(past peculiarity). There is no future singularity. For 
t — > +oo, the radius and the density tend to +oo and 
the pressure to — oo (Little Rip). The density increases 
exponentially rapidly, 
(iii) For n > — 2, 



R(t) 



2/(2- 



p* 



2-\n\ 



>-Kt 



2|n|/(2-|ra|) 



(A26) 



(A27) 



The universe starts at t = with a finite radius R = 1 
and a vanishing density (past peculiarity). At t = 0, the 
pressure vanishes for — 2<n<— 1, is finite for n = — 1 
and tends to — oo for n > — 1. In this last case, there is a 
past singularity of type II. On the other hand, there is no 
future singularity: For t — > +oo, the radius and the den- 
sity tend to +oo and the pressure to — oo (Little Rip). 
Actually, we can extend the solution to t < (except, 
maybe, in the case n > — 1 where the pressure diverges 
at t = 0). This describes a phase of contraction of the 
universe, corresponding to the solution of Eq. (|A2ip with 
e = — 1. Therefore, we obtain a model of bouncing phan- 
tom universe that collapses for t < (with decreasing 
density), disappears at t = (the density vanishes), and 
expands for t > (with increasing density). 

The evolution of the scale factor in these different cases 
is represented in Fig. IT3l 



Appendix B: Summary of all the possible cases 

The study of the polytropic equation of state ([T]) in 
cosmology is very rich. This is also the case for the study 
of polytropic distributions in stellar structure (3l| and 
in other areas of physics and biology j32l - l3^ . In this 
Appendix, we summarize all the results obtained in our 
series of papers and analyze the different singularities in 
terms of the classification of Gil . 



1. The case -1 < a < 1 

• In papers I and II, we have studied the case w > — 1. 
We have obtained the following results: 



(i) For n > and k < 0, the universe undergoes an 
early inflation. It starts from t = — oo with a vanishing 
radius a = and a finite density p*. Its radius increases 
indefinitely in time while its density decreases. There is 
no singularity. 

(ii) For n > and k > 0, the universe exhibits a past 
singularity of type III. It starts at t = with a finite 
radius a* and an infinite density p = +oo. For t > 0, 
its radius increases indefinitely in time while its density 
decreases. There is no future singularity. 

(iii) For n < and k < 0, the universe exhibits a Big 
Bang singularity. It starts at t = with a vanishing ra- 
dius a = and an infinite density p = +oo. The universe 
also undergoes a late inflation. Its radius increases to 
+00 as t — > +oo while its density decreases to a finite 
value p*. There is no future singularity. 

(iv) For ?i < and k > 0, the universe exhibits a Big 
Bang singularity. It starts at t = with a vanishing ra- 
dius a = and an infinite density p = +oo. The universe 
also exhibits a future peculiarity. Its radius increases to 
a finite value a* while its density decreases to zero p = 
(the universe "disappears"). For n < — 2, this peculiarity 
is reached in infinite time. For n > — 2, this peculiarity 
is reaches in a finite time t* (for n < — 1, there is a fu- 
ture singularity of type II because the pressure diverges 
at t = £*). For f* < t < 2t», the radius decreases to 
zero while the density increases to +oo. This leads to a 
Big Crunch singularity. These phases of expansion and 
contraction continue periodically (cyclic universe). 

• In this paper, we have studied the case w < — 1 
(requiring k < 0) corresponding to a phantom universe. 
We have obtained the following results: 

(i) For n > 0, the universe undergoes an early inflation. 
It starts from t = — oo with a vanishing radius a = and 
a finite density p». The universe also undergoes a future 
singularity of type III: At a finite time t* , its radius tends 
to a finite value a* while its density diverges p — > +oo. 

(ii) For n < 0, the universe exhibits a past peculiarity. 
Its radius starts from a finite value a* while its density 
vanishes p = 0. For n < — 2, this peculiarity occurs in 
the infinite past. For n > —2, this peculiarity occurs at 
t = (for n > — 1, there is a past singularity of type II 
because the pressure diverges at t = 0). The universe also 
undergoes a late inflation. Its radius increases to +oo as 
t — > +oo while its density increases to a finite value p*. 
There is no future singularity. Actually, the solutions 
with n > — 2 can be continued symmetrically for t < 
leading to models of bouncing phantom universe. 



2. The case a = — 1 

• In papers I and II we have studied the case w > 
— 1 (requiring k > 0). We have obtained the following 
results: 

(i) For n > 0, the universe exhibits a past singularity of 
type III. It starts at t — with a finite radius a* and an 
infinite density p = +oo. For t > 0, its radius increases 
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indefinitely in time while its density decreases. There is 
no future singularity. 

(ii) For n < — 2, the universe exhibits a Big Bang singu- 
larity. It starts at t = with a vanishing radius a = and 
an infinite density p = +oo. The universe also exhibits a 
future peculiarity. Its radius increases to a finite value a* 
while its density decreases to zero p = 0. This peculiarity 
is reached algebraically rapidly in infinite time. 

(iii) For n = —2, the universe starts from t = — oo with 
a vanishing radius a = and an infinite density p = +oo. 
The universe exhibits a future peculiarity. Its radius in- 
creases to a finite value a* while its density decreases 
to zero p = 0. This peculiarity is reached exponentially 
rapidly in infinite time. 

(iv) For — 2 < n < 0, the universe starts from t = — oo 
with a vanishing radius a = and an infinite density 
p = +oo. The universe exhibits a future peculiarity. Its 
radius increases to a finite value a* while its density de- 
creases to zero p = (the universe "disappears"). This 
peculiarity is reached in a finite time i* (for — 1 < n < 0, 
there is a future singularity of type II because the pres- 
sure diverges at t = £*). For t>t st , the radius decreases 
to zero while the density increases to +oo. This takes 
place in infinite time. 

• In this paper, we have studied the case w < — 1 
(requiring k < 0) corresponding to a phantom universe. 



We have obtained the following results (see also (l9l.[2]|): 

(i) For n > the universe starts from t = — oo with 
a vanishing radius a = and a vanishing density p = 0. 
This corresponds to a generalized past peculiarity. The 
universe also exhibits a future singularity of type III: At 
a finite time t* , its radius tends to a finite value a* while 
its density diverges p — > +oo. 

(ii) For n < —2, the universe exhibits a past pecu- 
liarity. It starts from t = — oo with a finite radius a* 
and a vanishing density p = 0. It also exhibits a future 
singularity of type I (Big Rip): At a finite time i*, its ra- 
dius and density are infinite. This singularity is reached 
algebraically rapidly. 

(iii) For n = — 2. the universe exhibits a past pecu- 
liarity. It starts from t = — oo with a finite radius a* 
and a vanishing density p = 0. Then, the radius and the 
density increase indefinitely (Little Rip). 

(iv) For —2 < n < 0, the universe exhibits a past 
peculiarity: It starts at t = with a finite radius a* and 
a vanishing density p = (for — 1 < n < 0, there is a 
past singularity of type II because the pressure diverges 
at t = 0). Then, the radius and the density increase 
indefinitely (Little Rip). Actually, the solution can be 
continued symmetrically for t < 0. This leads to a model 
of bouncing phantom universe. 
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